Weighted Trade Network in a Model of Preferential Bipartite Transactions 
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Using a model of wealth distribution where traders are characterized by quenched random saving 
propensities and trade among themselves by bipartite transactions, we mimic the enhanced rates 
of trading of the rich by introducing the preferential selection rule using a pair of continuously 
tunable parameters. The bipartite trading defines a growing trade network of traders linked by their 
mutual trade relationships. With the preferential selection rule this network appears to be highly 
heterogeneous characterized by the scale-free nodal degree and the link weight distributions and 
presents signatures of non-trivial strength-degree correlations. With detailed numerical simulations 
and using finite-size scaling analysis we present evidence that the associated critical exponents are 
continuous functions of the tuning parameters. However the wealth distribution has been observed 
to follow the well-known Pareto law robustly for all positive values of the tuning parameters. 

PACS numbers: 89.65.Gh, 89.75.Hc, 89.75.Fb, 05.70.Jk 



I INTRODUCTION 

In a trading society different traders trade among 
themselves. Thus the wealth distribution of the soci- 
ety dynamically evolves through this trading process. 
In the simplest possible situation pairs of traders make 
economic transactions. Such a mutual interaction can 
be looked upon establishing a connection between them. 
Consequently one can define a trade network where each 
trader is a node and a link is established between a pair 
of nodes when the corresponding traders make a mutual 
trading. In this paper we study the growth and the struc- 
tural properties of a trade network in the framework of 
a well known model of wealth distribution, namely the 
Kinetic Exchange Model (KEM) with quenched random 
saving propensities [l|, 0| . 

Over the last decade tremendous amount of research 
efforts have been devoted to study the structures, prop- 
erties and functions of different real-world as well as the- 
oretically defined model networks The key charac- 
teristic features of these networks include their small- 
world properties, which simply implies the existence of a 
very short global connectivity even when the sizes of the 
networks are extremely large [1]. Secondly, there are a 
large class of networks that are extremely heterogeneous. 
Their heterogeneity are quantified by their degree (num- 
ber of links meeting at a node) distributions. Usually 
such networks are observed to have power law decay of 
their degree distributions and are called Scale-free Net- 
works [5|. It has also been apparent that the links of 
many of these networks appear with a wide variation of 
strengths. In graph theoretic language the link strengths 
are called the 'weights' in general. Such weighted net- 
works have also been studied in the context of the pas- 
senger traffics of airport networks 0], international 
trade networks (ITN) etc. 

More than a century ago Pareto proposed that the dis- 
tribution of wealth x in a society to be P{x) ~ x~^'^~^'^'> 



[lo| . This form of distribution is generally known as 
Pareto distribution for a value oi v ^ 1 [ll|. Pareto 
suggested that j/ = 1 for the wealth distribution and it is 
known as the Pareto law. 

Application of the concepts of Statistical Physics to the 
wealth/income distribution in a society goes back to 1931 
when Saha and Srivastava had suggested that the form of 
the wealth distribution may be similar to the Maxwell- 
Boltzmann distribution of molecular speeds in an ideal 
Over the last few years renewed attempts 
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have been made using Statistical Physics methods. The 
main objective is to reproduce the recently collected in- 
dividual income tax data in different countries reflecting 
the wealth distributions. It has been observed that these 
data fit well to exponentially decaying functions for small 
wealths which however ends with power law tails in the 
large wealth regime. 

Dragulescu and Yakovenko (DY) [l3| modeled a bipar- 
tite trading between two traders using the analogy of a 
pair of gas molecules interacting through an energy con- 
served elastic collision in an ideal gas. Starting from an 
arbitrary distribution of individual wealths the system 
evolves through a series of bipartite trades to arrive at 
a stationary state where the wealth distribution assumes 
its time independent form. While the DY model produces 
only an exponentially decaying wealth distribution, later 
modifications were suggested for its improvement. This 
class of models are now referred as the KEMs lli,l|. 

In a KEM the society is considered as a collection 
of N traders, each having a certain amount of money 
equivalent to his wealth Xi,{i = 1,N} which he uses 
for mutual trades with other traders. Generally all 
traders are initially given an equal amount of money 
P{xi,t = 0) = S{xi — a). The sequential time t is the 
number of bipartite trades. A trade consists of two parts, 
(i) a rule for the selection of a distinct pair of traders i 
and j, {i =/= j) and (ii) a distribution rule for the ran- 
dom shuffling of their total money Xi + Xj between them. 
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Different KEMs differ from one another either in the se- 
lection rule or in the distribution rule or in both. In 
the DY model a pair of traders is selected with uniform 
probability. Their total money is then randomly reshuf- 
fled between them. In the stationary state the wealth 
distribution is P(a;) = exp{—x/{x))/{x) where the mean 



wealth (x) = a is usually set at unity Quite nat- 

urally a trader invests a only a part of his money in a 
trade and not his entire wealth. To incorporate this fact 
Chakraborti and Chakrabarti (CC) introduced a saving 
propensity A |15j] same for all traders. As a result the 
stationary state wealth distribution gets modified to a 
distribution with a single maximum which approximately 
fits to a Gamma distribution [IQ] . 



In a second modification of the DY model Chatterjee, 
Chakrabarti and Manna (CCM) assigned a quenched dis- 
tribution of saving propensities {Ajj so that each trader is 
characterized by his own Ai value [17| . Using the same se- 
lection rule as in DY model, the total money invested by 
a pair of traders after saving has been randomly shared 
between them. The system reaches a stationary state 
here as well but it sensitively depends on the precise 
values of A^s. The wealth distribution in the station- 
ary state after averaging over different realizations of the 
quenched disorder {Xi} yields a power law decay with a 
value of the Pareto exponent v ^ 1 [l^]. However, sub- 
sequent detailed analyses have revealed that the CCM 
model has many interesting features jl8| . For example, 
in contrast to the DY and CC models, CCM model is 
not ergodic. Therefore the wealth distribution is not 
self-averaging and the single trader wealth distribution 
is totally different from the over all wealth distribution 
of the whole society. Consequently the individual saving 
propensity factor A^ plays the role of an identification la- 
bel that determines the economic strength of a member 
in the society [3]. In fact, the wealth of a trader fluctu- 
ates around a mean value which depends very sensitively 
on the precise value of A^. Larger the value of A higher 
is the mean wealth. Truly the wealth distribution aver- 
aged over many uncorrelated quenched {A^} sets is the 
convolution of the individual members' wealth distribu- 
tions [3] ■ This overall distribution for the whole system 
exhibits Pareto law but not the individual member dis- 
tributions. The exponent v = 1 has been found to be 
exactly equal to unity in 1^ 2o| . 



In section II we describe our modification of the CCM 
model using the preferential selection rule. The wealth 
distribution of the modified CCM model has been de- 
scribed in section III. In section IV we define the asso- 
ciated trade network in terms of its nodal degrees and 
link weights. Section V presents the degree distribution 
and the weight distribution is discussed in section VI. We 
summarize in section VII. 



II THE MODEL 

It is a general observation that in a society the rich 
traders invest much more in trade and therefore take 
part in the trading process more frequently than the poor 
ones. To incorporate this fact in the CCM model that 
rich traders are preferentially selected more frequently 
with higher probabilities we introduce two parameters a 
and P in general, both > 0, to tune the preference differ- 
ent traders receive for their selections. We assume that 
the probability of selection of a trader is directly propor- 
tional to the a and /3-th power of its wealth. Therefore 
a pair of traders i and j {i ^ j) with money Xi and Xj 
are selected with probabilities 



TTi{t) (X Xi{t)" 



and 



7T,{t)^X,{tf. (1) 



When a = P = we get the ordinary CCM. When 
they are non-zero the rich traders are selected with larger 
probabilities. More rich a trader, higher is the probabil- 
ity that it will be selected for trading. Once a pair of 
traders i and j is selected, they save Ai and Xj fractions 
of their money and invest the rest amounts to the mutual 
trade. The total invested amount by both the traders is 
therefore 



6ij(t) = XiXi(t) + XjXjit), 



(2) 



where A = 1 — A. This amount is then randomly divided 
into two parts and received by them randomly: 



x,{t + l) ^ X,x^{t) + e{t)5,j{t) 

Xj{t + 1) = XjXj{t) + €{t)Sij{t). 



(3) 



where e(t) is a freshly generated random fraction and 
e = 1 - e. 

It is essential that all measurements are done once the 
system attains the stationary state. For that it is neces- 
sary that a number of transactions take place between ev- 
ery pair of traders, only then the mean wealths of traders 
attain their stationary values and fluctuate around them. 
Eqn. (1) states that for any {a,/3) > the richest and 
the next rich are the most probable pair and the poorest 
and the next poor are the least probable pair for selec- 
tion. Assuming the maximum wealth Xmax ^ N (with 
(x)=l) and the minimal wealth Xmin ^ 1/iV the relax- 
ation time can be estimated which is the typical time 
required for the poorest pair to make a trade. The poor- 
est is selected with a probability x^^^/Eixf . Approxi- 
mating the denominator by its maximum value we get 
{xmin/xmax)'^ ~ N"'^". Similarly the probability for 
the next poorest is N^"^^ and for the poorest pair as 
jY-2(a+/3) Therefore the time required for a trade be- 
tween the poorest pair T2 ~ N'^i°'+f^) (see below) and 
the relaxation time is several multiples of T2. Thus for 
any {a,P) > the relaxation time grows very rapidly 
with N. 
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FIG. 1: (Color online) (a) Plots of (a;(A))(l - A) vs. A for 
a = /? = (black), 1/2 (red), 1 (green), 3/2 (blue) and 2 
(magenta) for A'^ — 1204 (a increases from top to bottom), 
(b) The product function {x{X)){l - A)A"'''"' is plotted with 
A using x(a) = 0.15, 0.35, 0.57 and 0.80 for a = 1/2, 1, 3/2 
and 2 respectively using the same colors as in (a). 



At the early stage rich traders at the top level quickly 
take part in the trading but gradually the inclusion of 
relatively poor traders becomes increasingly slower. As 
a result the number of distinct traders taking part in 
the trading process grows very slowly. Effectively this 
implies that the system passes through a very slow tran- 
sient phase which is practically time independent. We 
call this state as the "quasi-stationary state (QSS)" . It is 
to be noted that in the following sections we present our 
numerical results for large system sizes in the QSS only. 
To ensure that the system has indeed reached the QSS in 
our simulations we keep track of the quantity S^x^ and 
collect the data only after no appreciable change of its 
mean value is noticed. We mostly analyse the symmetric 
a = f3 cases except for a few asymmetric cases. 



Ill WEALTH DISTRIBUTION 

For CCM it was observed that the average money 
of a trader (cc(A)) with saving propensity A diverges as 
A — » 1 [3]. Later it was shown that the divergence is 



ii, [23. This is simply be- 



like {xi){l — Xi) = constant 
cause had there been a trader with A = 1, he would not 
invest any money at all but always receives a share of 
the investments of the other traders! As a result this 
trader will eventually grab all the money of the society 
and, and this situation is similar to the phenomenon of 
condensation. 




FIG. 2: (Color onhne) (a) Wealth distribution P{x, N) vs. x 
for Q = /3 = 0, 1/2, 1 and 2 and for iV = 256 (black), 1024 
(red) and 4096 (blue) (A*' increases from left to right). The 
slopes of these curves yield v = 1.00(3) consistent with the 
Pareto law. (b) P{x, N) for {a, jS) = (00, 0) (black) and (0, 00) 
(red) for A'^= 1024 which almost overlapped, (c) P{x, N) for 
(a,/3) — (cx3,oo), the distribution is uniform followed by a 
hump due to transactions between the richest and the next 
richest traders only. 



In Fig. 1(a) we plot the quantity (a:(A))(l — A) vs. A 
for five different values of a = [3 = 0, 1/2, 1, 3/2 and 
2. For a = (3 = we see the horizontal line as observed 
in [21]. However for other a, f3 values the variations of 
the same quantity are far from being uniform and are 
monotonically increasing with A, their growth becoming 
increasingly faster with a. Therefore we try multiplying 
this function by A~^^"^ where x(c*) is a function of the 
parameter a. In Fig. 1(b) we plot (a:(A))(l — A)A^^^"^ 
vs. A using the same data of Fig. 1(a) using — 
0.15, 0.35, 0.57 and 0.80 for a = /3 = 1/2, 1, 3/2 and 2 
respectively. We get nearly uniform variations between 
A = 0.3 and 1. We assume that 



(x(A))(l - A)A~'^(") = constant. 



(4) 



If the distribution of A values is denoted by g{X) = 
constant, and since the wealth x and the saving propen- 
sity A are the two mutually dependent variables associ- 
ated with the same trader, their probability distributions 
follow the relation pjj] 



P(x)dx = g{X)dX. 



(5) 
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FIG. 3: The phase diagram in the positive quadrant of the 
(a, /3) plane. The Pareto law is valid in the entire region. The 
origin corresponds to the CCM model where the trade net- 
work is a random graph. At the corners {00, 0) and (0, 00) the 
richest trader trade in every transaction, so that the network 
is a star-like graph but the wealth distribution still follows 
Pareto law as shown in Fig. 2(b). However at the corner 
(00,00) the trade takes place between the top two richest 
traders and the network shrinks to a single dimer. The wealth 
distribution here is uniform followed by a hump as shown in 
Fig. 2(c). 



Differentiating Eqn. (4) with A one can find out the 
derivative dX/dx and substituting in Eqn. (5) one gets 



P(x) 



C 



[X-^ + (1 - A)xA 



-x-ii-i 



(6) 



For this equation we see that for large A the term within 
[. .] is of the order of unity. Therefore in this range P{x) ^ 
x^'^ as in the Pareto law. This is an indication that even 
for (a, 0) > 0, Pareto law holds good and in the following 
we present numerical evidence in support of that. 

The system is prepared by assigning uniformly dis- 
tributed random fractions for the saving propensities Ai 
to all N traders. Here A^s are quenched variables and 
therefore they remain fixed during the subsequent time 
evolution of the trading system. Consequently all ob- 
servable that we measured are averaged over different 
uncorrelated sets of the {Xi} values. While assigning the 
A values we first draw N uniformly distributed random 
fractions, but then to avoid the situation when Xmax is 
very close to unity by chance we scale them proportion- 
ately so that Xmax = 1 — 1/iV in every {A^} set. First a 
pair of values for (a, /?) is selected. Two types of initial 
wealth distributions are used: (i) Xi = I for all i and (ii) 
XiS are uniformly distributed random numbers with (x) 
— constant. The sequence of bipartite trading begins by 
randomly selecting pairs of traders using Eqn. 1. Once 
a pair is selected, their total individual invested amount 
Sij is calculated using Eqn. 2 and this amount is shared 
again between them using Eqn. 3. This constitute a sin- 
gle bipartite trading and the dynamics is followed over a 
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FIG. 4: (Color online) Growth of the trade network, (a) Plot 
of the average size of the giant component (sm(p, A'^)) with the 
link density p, for a = /3 = 1 and for N = 128 (black), 256 
(red), 512 (blue) and 1024 (pink), (A'^ increases from right 
to left). The inset shows a data collapse of the same plots 
with pN'^ , and 6 = 0.88. (b) The percolation link density 
Pc{N) is plotted with where 6 = 0.88, 0.92 and 1 for 

a — f3 = 1/2,3/4 and 1 respectively. The inset plots 6{a) 
with a. 



large number of such trading events. 

The wealth distribution changes with time from the 
initial distribution to more and more flat distribution. 
After a certain time the system passes through the 
quasi stationary state when no appreciable change in the 
wealth distribution is observed. It is also observed that 
the distribution is robust with respect to the precise val- 
ues of the parameters a and f3 used. In Fig. 2(a) the 
wealth distribution P{x, N) has been plotted with x for 
four sets of parameter values namely, a = /? = 0, 1/2, 1 
and 3/2 and for three system sizes N = 256, 1024 and 
4096. Apart from slight fluctuations the four curves for 
a given system size nearly overlap on one another. On 
a double logarithmic scale the slopes of the curves give 
an average estimate for ly = 1.00(3) consistent with the 
Pareto law as observed in the CCM [l3|. This indicates 
that the wealth distribution is robust with respect to the 
parameter values in this region. The non-zero values of a 
and P only controls the frequencies with which different 
traders are called for trading. 

Next we consider the case when one of the two pa- 
rameters {a, (3) is infinity and the other one is zero. If 
a = 00 the richest trader is always selected as the first 
trader. The other trader is selected among the other 
A'^ — 1 traders with uniform probability. As shown in Fig. 



5 




I I I I I , I I I I I I I , I I I I I , I 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

a 



FIG. 5: (Color online) (a) The finite-size scaling of the degree 
distributions P{k,N) for a = /? = 1, for iV = 256 (black), 
1024 (red) and 4096 (blue), and for (fc) = 1. Direct mea- 
surement of slopes give 7^(1) = 2.18(3). The best data col- 
lapse corresponds to 7?fc(l) = 1.62 and Cfc(l) = 0.75 giving 
7fe(l) = 2.16(3). (b) The plot of 7fc(Q) vs. a. 

2(b) here also we see that the Pareto law holds good. For 
a = oo and for finite /3 first the richest trader is selected 
and then the second trader is selected with probability 
oc . We observe numerically that here also Pareto law 
works very well. 

However the situation is very different when both 
(a,/9) take very large values. In this situation almost 
always only the rich traders are called for transactions. 
The system passes through an extremely long QSS and 
the number of traders taking part in trade does not in- 
crease at all. For example in the limiting case of (a, (3) 
— (cxD, cxd) it implies that always only the richest and 
the next richest traders are selected for transactions with 
probability one but not any other trader. If their wealths 
are very high then the trading will be limited only be- 
tween them. Therefore the wealth distribution for the 
single Aj set has two very high peaks and wealths of all 
other traders are small and uniformly distributed. Con- 
sequently the quench averaged wealth distribution is uni- 
form throughout followed by a hump at the highest value 
of wealth (Fig. 2(c)). A systematic analysis with many 
different (a, /3) pairs leads us to conclude that Pareto law 
holds good in the positive quadrant of the entire (a, (3) 
plane. 

In Fig. 3 we exhibit this behavior in the positive quad- 
rant of the (a, P) plane where Pareto law is valid and the 
limiting points are marked by circles with their character- 
istics. The origin at (a = 0, /3 = 0) represents the CCM 
model where traders are selected randomly with uniform 
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FIG. 6: (Color online) Probability distribution P(w, a) of the 
link weights, (a) Plots for the A'^-clique graphs with =64 
and for a— (black), 1/4 (green), 1/2 (blue), 3/4 (magenta) 
and 1 (red) and (3 =1 always, (b) Plots for (fc) — 5 with 
a — P — 1 and for the system sizes 128 (black), 256 (red) and 
512 (blue). Direct measurement of slopes gives 2.52, 2.53 and 
2.51 respectively. 

probabilities. As explained below, the trade network cor- 
responding to this point is a random graph (RG). As ex- 
plained before that at the two corners (oo, 0) and (0, oo) 
the richest trader always participates in every transac- 
tion. Therefore the corresponding trade networks have 
star-like structures. In the last corner of (oo, oo) the 
trade takes place only between the richest and the next 
rich traders and therefore the graph reduced to a single 
dimer only. 

IV THE TRADE NETWORK 

One can associate a network with this trading system. 
Each trader is a node of the network. Initially the net- 
work has only N nodes but no links. First the system 
is allowed to reach the QSS and then the network starts 
growing. Every time a pair of traders makes a trade for 
the first time, a link is introduced between their nodes. 
There after no further link is added between them irre- 
spective of their subsequent trades and they remain con- 
nected with a single link. As the system evolves more and 
more new traders take part in the trading dynamics and 
consequently the number of links grow in the network. 
For a — P — the growth of the network is exactly the 
same as that of the random graph, however it is much 
different when {a, /?) > 0. Since the rich nodes are pref- 
erentially selected they acquire links at a faster rate than 
the poor nodes. The degree ki of the node i is the num- 
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FIG. 7; (Color online) (a) Nodal strength distributions 
P{s,N) with strength s with a = /3 = 1, for iV = 256, 1024 
and 4096 and for (fc) — 5. Direct measurement of the slopes 
of these curves gives 7s — 2.50(5). (b) The Finite-size scal- 
ing analysis of the data in (a) gives 77s = 1.64 and (s ~ 0.67 
estimating the value of 7s = »7s/Cs = 2.45(5). 
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FIG. 8: (Color online) (a) Plot of nodal strength {s{k)) with 
degree fc for = 2" and for a = /3 = 1/2 (black), 3/4 
(green), 1 (blue), 5/4 (magenta) and 3/2 (red) (from bottom 
to top). The slopes of these plots estimate </>(«) shown in 
the inset, (b) Plot of average wealth {x{k)) with degree k 
for iV = 2" and for a = /3 = 1/2 (black), 3/4 (green), 1 
(blue), 5/4 (magenta) and 3/2 (red) (from bottom to top). 
The slopes of these plots estimate /i(a) shown in the inset. 



ber of distinct traders with whom the i-th trader has ever 
traded. The dynamics is continued for a certain time T 
till the average degree (k) of a node reaches a specific 
pre-assigned value. 

In general there are two characteristic time scales in- 
volved. At the early stage the network grows with mul- 
tiple components with different sizes. At time Ti the 
growing network becomes a single component connected 
graph. A second time scale is T2 when the whole network 
is a A'^-clique graph in which each node is linked to all 
others, which means each trader has traded at least once 
with all others. Unlike random graphs the growth of the 
network is highly heterogeneous and the rich traders have 
much larger degrees than the poor traders. Since poor 
traders are selected with low probabilities they take much 
longer times to be a part of the network. Consequently Ti 
is found to be very large and of the same order as T2. Nu- 
merically it is easier to calculate T2, one keeps track of the 
number of distinct links and stops only when this number 
becomes just equal to [N{N — l)]/2. On the other hand 
to calculate Ti one follows the growth of the giant com- 
ponent and stops when the giant component covers all N 
nodes. A Hoshen-Kopelman cluster counting algorithm 
[2^ is used to estimate the size of the giant component. 
For the ordinary CCM with a = f3 = since both traders 
are chosen with uniform probability, the generated graph 
is a simple Erdos-Renyi random graph characterized by 
a Poissonian degree distribution [23] . 



The growth of the giant component is studied with 
increasing number of links n in the in the trade network. 
The average fraction of nodes in the giant component 
is denoted by {sm{p,N)) which is the order parameter 
in this percolation problem. This has been plotted in 
Fig. 4(a) using a semi-log scale with link density p — 
n/[N{N — l)]/2 in the network. Four curves shown in 
this figure correspond to iV = 128, 256, 512 and 1024 for 
a = f3 = 1, the system size increasing from right to left. 
The inset shows that a data collapse can be obtained by 
scaling the p axis by a factor with 9 = 0.88. The 
critical density of percolation transition Pc{N) is defined 
as that particular value of p for which the average size of 
the giant component {sm{p,N)) = 1/2. In Fig. 4(b) we 
show that how the critical percolation threshold pc{N) 
depends on by plotting it with for a = /3 = 1/2, 
3/4 and 1. It has been observed that the exponent 9 (a) 
is dependent on a in general and in the inset of this figure 
we plot 9{a) vs. a. We see that for a < 1/2, x{a) = 1 but 
for Of > 1/2, x{a) gradually decreases. For Erdos-Renyi 
random graphs it is known that a: = 1 and therefore this 
result gives an indication that the trade network seems 
to be different from random graphs for a = [3 > 1/2. 
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V DEGREE DISTRIBUTION 

The degree distribution has been studied similar to 

random graphs. Wo keep track of the average degree (k) 
of the network which is related to the number of hnks n 
of the network as n = {k)^. First the degree distribution 
has been studied for (fc) = 1 and for different system sizes. 
For an assigned set of values of (a,/?), for a given set of 
values for the saving propensities and for a specific 
value of N once (fc) = 1 is reached we calciilate the degree 
distribution considering all components of the network 
on the same footing. The network is then refreshed by 
removing all links and a second network is constructed 
and so on. The dynamics is continued for the same values 
of the parameters and the same set of till a large 
number of networks are generated and their mean degree 
distribution is calculated. The entire dynamical process 
is then repeated with another uncorrelated set of {Aj}s 
and the degree distribution has been averaged over many 
such sets. 

In Fig. 5(a) we show the finite-size scaling plot of 
the average degree distribution P(fc, N) vs. k for a = 
/3 = 1 and for A^^ = 256, 1024 and 4096. On a double 
logarithmic scale all three curves show quite long scaling 
regions followed by humps before the cut-off sizes of the 
degree distributions. The cut-off sizes of the distributions 
shifts to the larger values of k approximately by equal 
amounts on the double-log scale when the system size 
has been enhanced by the same factor. From the direct 
measurement of slope in the scaling region we estimate 
7(1) = 2.18(3). Almost the entire degree distribution 
obeys nicely the usual finite-size scaling analysis and an 
excellent collapse of the data is observed confirming the 
validity of the following scaling form: 

P{k, N) oc iV-''(«)a[A;/A^^<"'] (7) 

where the scaling function G{y) has its usual forms like, 
S{y) ~ y~'>'(") as y and G{y) approaches zero 
very fast for y >> 1. This is satisfied only when 
7(0;) = ({a)/ri{a). The exponents rj{a) and Ci^^) fully 
characterize the scaling of P{k, N) in this case. To check 
the validity of the equation we attempted a data collapse 
by plotting P{k, N)N'^'^^'^ vs. fc/7V^(i) by timing the val- 
ues of 77(1) and C{1). The values obtained for best data 
collapse are ?7fc(l) = 1.62 and Cfc(l) = 0.75 implying that 
in the infinite size limit P(fc, 00) ^ k~'^^^"> with 7(1) = 
2.16(3). Tuning a and (3 to other values it is observed 
that the degree distribution exponent 7 does depend on 
these two parameters. In Fig. 5(b) we show a plot of 
Jk{ct) with a which decreases to « 2 at a = 2. 

VI THE WEIGHTED NETWORK 

Within a certain time T a large number of bipartite 
trades take place between any arbitrary pair of traders. 



The total sum of the amounts Sij invested in all trades 
between the traders i and j in time T is defined as the 
total volume of trade Wij = T,'^6ij. Therefore wtj is re- 
garded as the weight of the link (ij). The magnitudes of 
weights associated with the links of the trade network are 
again found to be highly heterogeneous. This is primar- 
ily because within a certain time T a rich pair of traders 
trade many more times than a rich-poor or a poor-poor 
pair. In addition the invested amounts depend on the 
mean wealths (xi) of the traders involved as well as their 
saving propensity factors Aj. The probability distribu- 
tion P(w, N) of the link weights are calculated when the 
average degree (k) reaches a specific pre-assigned value. 
As before, this distribution has also been averaged over 
many weighted networks for one {Aj} set and then fur- 
ther averaged over many imcorrelated {Xi} sets. 

First we studied the case when the trade networks is 
a A/'-clique graph, i.e., when each trader has traded with 
all other traders at least once. Here each node has same 
degree i.e., P{k) = 6{k - {N - 1)) and {k) = N - 1. The 
required time T2 increases rapidly with N as described in 
section II and we could study small system size A^ = 64 
only. The distribution has a very long tail and therefore 
we used a lin-log scale for plotting. In Fig. 6(a) we show 
the plots of P{w, a) with ln(w) for different values of a = 
0, 1/4, 1/2, 3/4 and 1 and /3 = 1. Each curve is asymmet- 
ric and has a single maximum. The position of the peak 
shifts towards larger values of ln(w) as a increases. If Fig. 
6(b) a similar plot has been shown for (k) = 5 for three 
network sizes N = 128, 256 and 512 and for a = /? = 1. 
On a double-logarithmic scale each curve has a consider- 
able straight portion. This indicates a power law decay 
like P{w,N) oc w"'''™. The corresponding slopes give 
estimates for the exponent 7^ as 2.52, 2.53 and 2.51 for 
the three system sizes respectively, so that on the average 
7^ = 2.52(3). 

The strength of a node Si = 'SjWij where j runs over 
all neighbors ki of i, is a measure of the total volume of 
trade handled by the i-th node. Nodal strengths varies 
over different nodes over a wide range. We first study the 
probability distribution of nodal strengths. In Fig. 7(a) 
the strength distribution P{s, N) has been plotted for the 
average degree (k) =5, for a = /3 = 1 and for the network 
sizes 256, 1024 and 4096. Extended scaling regions at the 
intermediate regions of the curves indicate that P{s, N) 
also follows a power-law decay function P(s, N) ~ s~^'' 
in the limit of N —>■ 00. Direct measurements gives an 
estimate of 7s (1) ~ 2.5. In Fig. 7(b) we try a similar 
finite size scaling of the same data giving r/s(l) = 1.64 
and C(l) = 0.67 giving 7^(1) = 2.45(5). 

Quite often weighted networks have non-linear 
strength-degree relations reflecting the presence of non- 
trivial correlations, example of such networks are the air- 
port networks and the international trade network. For a 
network where the link weights are randomly distributed, 
the {s{k)) grows linearly with k. However a non-linear 
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growth like (s(fc)) ~ k'^ with (j> > 1, exhibits the pres- 
ence of non-trivial correlations. In Fig. 8(a) we plot the 
variation of {s{k)) vs. k for a system size A'' = 16384 and 
for different values ofa = /3 = l/2 (black), 3/4 (green), 
1 (blue), 5/4 (magenta) and 3/2 (red) (from bottom to 
top) . The slopes of these plots give estimates for the ex- 
ponent (/'(a) which gradually increased with a and the 
variation has been plotted in the inset. In the same con- 
text we also studied how the mean wealth of a trader de- 
pends on its degree. The mean wealth of a trader {x{k)) 
has been plotted in Fig. 8(b) with its degree k for the 
same system sizes as in Fig. 8(a) and for the same values 
of parameters. A power law growth has been observed 
for all values of a: {x{k)) ^ k^^"^ The slopes of these 
plots give estimates for the exponent /Li(a) which has been 
plotted in the inset of Fig. 8(b). 

IV SUMMARY 

To summarize we have studied the different structural 
properties of a trade network associated with the dy- 
namical evolution of a model of wealth distribution with 
quenched saving propensities. In this model distinguish- 
able traders make preferential bipartite trades among 
themselves and in this way create links. They are se- 
lected for trade preferentially using a pair of continuously 
tunable parameters, where the rich traders are picked up 
more frequently for trade than poor traders. This creates 
huge heterogeneity in the system which has been reflected 
in the power-law distributions of the nodal degree and the 
link weight distributions measuring the volumes of trade. 
We present numerical evidence that the associated indi- 
vidual wealth distribution follows the well known Pareto 
law robustly for all positive values of the parameters. 

We thankfully acknowledge P. K. Mohanty, A. Chat- 
terjee and B. K. Chakrabarti for discussion and critical 
reading of the manuscript. 
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